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R.J Tweedie 


1. INTRODUCTION 

The application of the RHMC algorithm [p to lattice QCD with Ginsparg-Wilson quarks 
together with the computational power of QCDOC has resulted in the ability to generate, and per¬ 
form measurements on, dynamical 2-1-1 flavour configurations for the first time. However, before 
embarking on a large-scale and (potentially) computationally costly production run, it is necessary 
to explore the available parameter space with (several) smaller and, therefore, computationally 
cheaper initial trials. This work focuses on the recent ensembles generated on the QCDOC ma¬ 
chines for exactly this purpose. 

2. SIMULATION PARAMETERS 

The analysis was carried out on 2-1-1 flavour domain wall fermion configurations generated on 
the QCDOC machines. The standard domain wall Dirac operator |§] and Pauli-Villars field with 
the action introduced in |^] was used. The gauge action was defined by 


Sg[U] = -^ (l-8ci) £ P[U],,^v + ci £ RpUy 


( 2 . 1 ) 



where P[U]x.^v and/?[f/]x,^v represent the real part of the trace of the path ordered product of links 
around the 1x1 plaquette and 1x2 rectangle, respectively, in the , V plane at the point x, and p = 
6/g^ with g the bare coupling constant. For the DBW2 gauge action 0, ||], the coefficient ci was 
chosen to be —1.4069, whereas for Iwasaki ci = —0.331. In the DBW2 case the ensembles 
have three different j3 values: 0.72, 0.764 and 0.78, while for the Iwasaki case they have two: 2.13 
and 2.2. All ensembles were generated using the RHMC algorithm |[T]] with a trajectory length of 
0.5, volume of 16^x32, fifth dimension length of 8, and a domain wall height of 1.8. Two different 
ensembles were generated at each j3 value, one with a light isodoublet with mass tUud = and 
one where the light sea quark masses are both equal to our rough estimate of the strange quark 
mass, nis- An additional DBW2 ensemble with a light isodoublet mass niud = \ms at j8=0.72 was 
generated. The number of trajectories in each ensemble and the number of configurations used in 
the analysis are shown in table [I[ The ensembles had aniud = 0.01, 0.02 or 0.04 and anis = 0.04. 

Measurements were made with up to four valence quark masses in the range 0.01 to 0.04 on 
each ensemble and, on some of the ensembles, correlators were measured with sources on multiple 
time planes to improve our statistics. Several types of smearing were used, in particular, point 
sources, wall sources, and hydrogen-like wavefunction smearing where one or two of the quark 
propagators in a meson correlator were smeared at the source. In general, simultaneous fits to local 
and smeared correlators were performed throughout. This analysis aggregates to more than thirty 
thousand trajectories and more than one hundred and twenty thousand measurements. 

The integrated auto-correlation time for the pseudoscalar meson correlator was calculated on 
several of these ensembles and found to be of order 100 trajectories. The correlators were over¬ 
sampled and averaged into bins of size between five and ten depending on whether the separation 
between measurements was five or ten trajectories. Our observed errors stabilised with bins of 
this size in good agreement with the calculated integrated auto-correlation time. A full correlated 
analysis was then performed with the binned data as input. 
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Action 

P 

mud 

Mr 

aruud 

Ntraj 

^meas 


DBW2 

0.72 

1.0 

0.04 

3395 

475 


DBW2 

0.72 

0.5 

0.02 

6000 

1000 

0.0106(1) 

DBW2 

0.72 

0.25 

0.01 

2340 

269 


DBW2 

0.764 

1.0 

0.04 

1615 

150 

0.00546(7) 

DBW2 

0.764 

0.5 

0.02 

1800 

88 

DBW2 

0.78 

1.0 

0.04 

1620 

113 

0.00437(6) 

DBW2 

0.78 

0.5 

0.02 

1505 

202 

Iwasaki 

2.13 

1.0 

0.04 

2380 

139 

0.0104(2) 

Iwasaki 

2.13 

0.5 

0.02 

2450 

146 

Iwasaki 

2.2 

1.0 

0.04 

4565 

713 

0.0065(1) 

Iwasaki 

2.2 

0.5 

0.02 

3175 

435 


Table 1: RHMC 2+1 flavour datasets and the value of the residual mass in the chiral limit. 

3. RESULTS 

3.1 The Residual Mass 

The residual mass is a measure of the violation of chiral symmetry [Q]. In our ensembles 
the length of the fifth dimension is relatively short, hence there is a significant left-right coupling 
between the quark fields on opposite walls. The residual mass was calculated from the ratio of the 
point-split pseudoscalar density J5 at the middle of the fifth dimension to the pseudoscalar density 
P built from the fields on the walls 


^x,y{J5{y,t)P{x,0)) 

A good signal for the value of am^es was observed, as can be seen in figure [I] (leff). 
Chiral extrapolations were performed by defining the quark mass 


(3.1) 


airiq = a{m f + f)), (3.2) 

where anif is the valence quark mass and is the residual mass measured using quark 

propagators with that quark mass, and taking aniq —> 0 using only the points where the valence 
quark mass is equal to the u and d quark masses in the sea. Other than in the DBW2 {3=0.72 
case, where a linear fit to three points was performed, straight lines were drawn through the two 
available points. The values of the residual mass obtained in the chiral limit are shown in table |^. 
These values correspond to a residual mass of ~5-30 MeV. 


3.2 Light hadron masses 

Pseudoscalar and vector meson masses are extracted by performing simultaneous double cosh 
fits to extract both the ground and first excited states. 

Chiral extrapolations were performed in an analogous way to the residual mass by extrapolat¬ 
ing the results to niq = 0. In the case of the pseudoscalar mass this takes the form 

{ampsf = B{amf=^‘^‘^)+A. (3.3) 
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Figure 1: LEFT: Residual mass for j5 =2.13 ensemble with = jm^. Different colours correspond to 
whether zero, one or two quark propagators are smeared at the source. RIGHT: Pseudoscalar meson decay 
constant using two different methods versus quark mass for the DBW2 j3=0.72 ensembles. 


Given the low statistics, an acceptable slight deviation of A from zero, typically less than 2a, is 
observed. The lattice spacing was obtained from nip in the chiral limit by performing a linear chiral 
extrapolation of the vector meson mass 

amv = +D. (3.4) 


This, together with the physical kaon mass and substituting —> amq^ + 

allowed the evaluation of the strange quark mass by setting aniq^ = amud = 0 and anig^ 


in eq.( [3.3| ), 
= am^. 


3.3 Pseudoscalar meson decay constants 


The pseudoscalar meson decay constant defined by 


afps 


Za{0\A4\P) 

amps 


(3.5) 


was calculated in two ways. In the first method was obtained from the axial Ward-Takahashi 
identity 

Za (dpApff) = 2{amq + {Pff) (3.6) 


which can be expressed in terms of the pseudoscalar meson correlator, Cpp, and the pseudoscalar 
meson axial correlator, Cpa^- ( 0 |A 4 |P) cancels when Za is substituted in eq.(^) and hence only 
Cpp and the value of the residual mass are required in order to evaluate fps 


afps 


2{am f + am,-gfi) {0\P\P) 
{ampsY 


(3.7) 


The red symbols in figure |T] (right) show the chiral extrapolation of fps for the DBW2 j3=0.72 
ensembles. 
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In the second method, we calculated the value of Za explicitly from a ratio designed to remove 
^(a) and suppress ^{a^) lattice artifacts ||^. 


C(f + i)+C(f —i) C(t + i) 

4l(0 +(L(0+L(t + 1)) 


(3.8) 


where C{t) and L{t) are the correlators of the pseudoscalar density with the partially conserved 
and local axial currents respectively. A simultaneous fit to both point-point and wall/smeared-point 
correlators is used to extract this ratio. The axial-axial correlator was then used in combination 
with the value of Za to calculate the pseudoscalar meson decay constant. The chiral extrapolation 
of fps using this method is shown by the blue symbols in figure |I] (right). We see good agreement 
between both methods. 

Using the lattice spacing and the strange quark mass calculated from the vector meson and 
pseudoscalar meson chiral extrapolations, we are able to calculate fx and extract the ratio 

JJt 

3.4 Scaling 

The ensembles generated for this analysis with several different lattice spacings and two dif¬ 
ferent gauge actions have ^(a^) discretisation errors. Figure ^ (left) shows the ratio ^ plotted 
against lattice spacing squared. The different colour symbols correspond to the different gauge ac¬ 
tions, Iwasaki in blue and DBW2 in red. We have set the lattice spacing in two different ways. The 
closed symbols correspond to setting the lattice spacing from the p meson mass, while the open 
symbols are from setting the lattice spacing from f^. Note that the errors are still quite large. The 
furthest two points to the right are the DBW2 p=0J2 points where we have been able to perform a 
chiral fit compared to drawing straight lines through the points so their errors are better estimated. 
Even with the large and crude error estimate for the other results, it may be concluded that there 
is an indication of scaling and, if anything, scaling may slightly better for the Iwasaki gauge ac¬ 
tion than for DBW2. Figure ^ (right) shows a similar plot to figure ^ (left), for ^ and similar 
behaviour is seen. 


4. SUMMARY 

Several ensembles have been generated on the QCDOC machines using the domain wall 
fermion formulation with two different gauge actions, several j3 values and multiple sea quark 
masses. These ensembles have a relatively small volume and limited statistics as they were primar¬ 
ily generated to search parameter space for a larger production run. A fifth dimension size of eight 
produces a residual mass larger than would be acceptable for such a production run. Even with 
these drawbacks it is still possible to calculate the light hadron spectrum and pseudoscalar decay 
constants obtaining results which are consistent with experiment and scale within large errors. At 
present we are running an ensemble with a larger volume of 24^x64 and a fifth dimension size of 
16 using the Iwasaki gauge action at j3 = 2.13. 
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Figure 2: Scaling behaviour of LEFT: ^ and RIGHT: for all ensembles. Different colours correspond 

to different gauge actions. Open and closed symbols correspond to setting the lattice spacing from or nip 
in the chiral limit respectively. 
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